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We show that both confined atoms and electron-atom scattering can be described by a unified 
basis set method. The central idea behind this method is to place the atom inside a hard potential 
sphere, enforced by a standard Slater type basis set multiplied by a cutoff factor. For confined atoms, 
where the wall is placed close to the atomic nucleus, we show how the energy of the highest occupied 
atomic orbital and the static polarizability of helium and neon atoms evolve with the confinement 
radius. To our knowledge, these are the first confined atom polarizability calculations that include 
correlation, through the use of time-dependent density-functional theory. By placing the atom in a 
large spherical box, with a wall outside the electron density, we obtain scattering phase shifts using 
a recently developed method [M. van Faassen, A. Wasserman, E. Engel, F. Zhang, and K. Burke, 
Phys. Rev. Lett. 99, 043005 (2007)]. We show that the basis set method gives identical results to 
previously obtained phase shifts for e-H and e-He + scattering. 



I. INTRODUCTION 

In this paper we study the behavior of atoms under 
the influence of a hard spherical potential wall. We posi- 
tion this wall both close to the atomic nucleus, inside the 
range of the free atom density, and far from the nucleus, 
outside the range of the free atom density. 

The first of these two cases describes a confined 
atom. The effects of spatial confinement are important 
in the description of several chemistry related phenom- 
ena. Some examples are the catalytic behavior of sur- 
faces and cages, such as zeolites, and the properties of 
electronic and excited states in quantum dots. For ex- 
amples of confined systems see Refs. p], and references 
therein. There are several methods available to impose 
the boundary condition of an infinite spherical potential 
on the system. One of these is to explicitly impose Dirich- 
let boundary conditions. This method was previously 
used to obtain low-energy electron-atom scattering phase 
shifts from (time-dependent) density-functional theory 

((TD)DFT) H H la]. We can also im P licitl y impose 
the boundary condition by requiring all basis functions 
to vanish at the wall boundary. In the classic work of 
Ludena 0, 0] a basis set for confined atoms was gener- 
ated by multiplying standard Gaussian and Slater type 
orbitals (GTO, STO) with a cutoff factor. More recently, 
other basis set methods were developed. One of these 
optimizes the orbital exponents by fitting them to cutoff 
orbitals Q and another generates a confining potential 
directly from the basis functions @. 

In the second case, where the wall is outside the free- 
atom density, the bound states of the system are hardly 
affected by the presence of the wall. Here the purpose 
of the wall is to discretize the continuum. We can ob- 
tain low-energy electron-scattering properties of atoms 
from these discretized continuum states 0, OH [H| • Low- 
energy electron-scattering has received a lot of inter- 
est since the discovery of Sanche and co-workers that 
low-energy electrons can cause single and double strand 
breaks in DNA molecules [l2j |. Such large systems are 



hard to describe by present scattering theories, although 
progress has been made in the last years pjl, [3, El- 
To provide an alternative method, an electron-molecule 
scattering method u sing (TD)DFT is currently under de- 
velopment 0, i, H3, Ha • Use of the TDDFT method in 
conjunction with discretized continuum states has proven 
to be an efficient method to obtain electron-atom scat- 
tering phase shifts 0, [Bj] . 



In this paper we use basis functions to describe the 
discretized continuum states. Finite basis sets are of- 
ten used to describe (discretized) continuum functions, 
even though they might not seem to be the most logical 
choice. The wish to use existing bound state methods 
and treat bound and continuum states on equal footing 
fuels the need for such basis set methods. Often, GTOs 
are employed [H, Gil US EH, [22[ , but the use of B-spline 
functions is also gaining popularity (see Ref. [23[ and ref- 
erences therein). We use the method of Ludena and mul- 
tiply STO basis functions with a cutoff factor to impose 
a hard wall outside the atom. We show that such con- 
fined basis functions are not only suitable for describing 
confined atoms, but also to obtain low-energy electron 
scattering properties. We hence obtain a unified basis 
set method for both confined atoms and electron-atom 
scattering. 



In the first part of this paper we introduce the ba- 
sis set for confined atoms and electron-atom scattering. 
Next, we explain how the properties studied in this pa- 
per, the energy of the highest occupied molecular orbital 
(HOMO) and polarizability for confined atoms and the 
electron scattering phase shifts, can be obtained from 
(TD)DFT. In the computational details section we ex- 
plain in detail how the confined STO basis functions can 
be made suitable for obtaining continuum states. Finally, 
we show results for confined He and Ne atoms, and elec- 
tron scattering from H and He + . 
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II. THEORY 

In this section we explain how we place a hard wall 
around an atom by introducing a cutoff function in the 
basis set. Next, we show how a particular choice of basis 
set parameters leads to a basis set suitable for obtain- 
ing discretized continuum functions. Finally, we explain 
how the polarizability and scattering phase shifts can be 
obtained from time-dependent density-functional theory 
(TDDFT). 



A. Boundary condition Slater type orbitals 

Most practical solutions of the Schrodinger equation 
involve the use of basis functions of some kind. To obtain 
the lowest bound states, generally a linear combination of 
atomic orbitals (LCAO) expansion, employing GTOs or 
STOs centered on the different atomic nuclei, is applied. 
GTOs are far more popular than STOs since the required 
integrals can be easily calculated analytically. A draw- 
back of GTO basis sets is that many GTO functions are 
necessary to describe the cusp behavior at the nucleus. 
In practice, contracted GTO basis sets are used, which 
reduce the computational effort that an uncontracted set 
would require. STOs have the advantage of possessing 
the required cusp behavior as well as the appropriate 
long-range decay. This allows the construction of high- 
quality basis sets with a relatively small number of func- 
tions. STO basis sets can be efficiently implemented in 
DFT codes when a density fitting procedure is used to ob- 
tain the Coulomb integrals [24| . A type of basis set that 
more recently gained popularity in atomic and molecu- 
lar calculations are basis splines, or B-splines. B-splines 
have obtained widespread use in atomic physics [23| be- 
cause they can be adapted readily to the problem under 
investigation. This means that a B-spline based code, op- 
timized to the study of bound states, can be adapted to 
continuum-state problems with minimal effort. Only the 
knot set, the small set of parameters defining the basis 
set, needs to be changed. In the next section we compare 
the continuum basis functions we develop on the basis of 
STOs to a typical B-spline basis. In this section we focus 
on the use of STO basis functions to describe confined 
atoms. 

A hard wall boundary is imposed by multiplying the 
basis functions with a cutoff function. The use of cutoff 
functions in conjunction with STOs has been first pro- 
posed by Ludeha [6, 7], who suggested the following form 
of the STOs, which we call boundary STOs (BSTOs) 



= Nf k (r;R)R n (r)Y lrn (0,< 
= NF BSm (r)Y lm (0^). 



(1) 



In this expression TV is a normalization factor and n, 
and m are the usual quantum numbers. The >^ m (#,(/>) 
are spherical harmonics. The radial part R n (r) is the 



same as for ordinary STOs 

Rn(r) = r n - x e~ aT , 



(2) 



where a is the orbital exponent. The cutoff function 
/fe(r; R) is defined by Ludeha as, 



fk(r;R) 







r <R 
r > R 



(3) 



This function becomes zero at a sphere of radius i?, fixing 
the boundary condition at R. For k > 2 the first deriva- 
tive of a BSTO is also zero at i?, and the function and its 
derivative are continuous across R. We choose k = 1, this 
choice has been found to lead to the most accurate results 
for confined atoms [9[. For k = 1 the first derivative has 
a jump discontinuity at i?, which obviously introduces 
problems when calculating first and second derivatives 
numerically. Therefore we take the analytic expressions 
of all necessary derivatives, avoiding numerical problems. 



B. A BSTO set for obtaining pseudo continuum 
states 

The BSTOs turn out not only to be suitable for calcu- 
lating properties of confined atoms, but also to calculate 
discretized continuum states and scattering phase shifts. 
Ordinary STO or GTOs are in general not suitable for 
obtaining scattering states because of their exponential 
decay, which is notably different from the oscillating be- 
havior of continuum states. At a first glance, the same 
seems to be true for BSTO basis functions, but they turn 
out to be quite suitable for the job at hand. The main 
reason for this difference is the use of the cutoff function. 
A normal STO (or GTO) will always become exponen- 
tially small at large radii and therefore is not able to 
accurately describe an oscillating continuum state. We 
can choose the exponential coefficients of our BSTO basis 
such that some of our basis functions are not exponen- 
tially small at the radius of the confining sphere R. The 
cutoff function then ensures that these basis functions 
are zero exactly at i?, but nonzero before. 

In this section we explain how we obtain the param- 
eters that define our BSTO basis. We show why this 
choice of parameters is especially suitable for scattering 
calculations and show their similarity to B-spline func- 
tions, which have already proven their merit in scattering 
calculations (23[. 

Our BSTO basis consists only of functions with n > 2. 
These functions have maxima located at some value of r, 



c . At r max we have 



dF 



BSTO 



(r) 



dr 



= °=R e 



-2 

' max 



'» [(n - i)i?<: 
(n + aR)r%£ + a< ax ] . (4) 



From this equation we obtain a relation between the un- 
known parameters n and a and the known parameters R 
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FIG. 1: Full lines: BSTO basis consisting of 5 functions 
with ristart = 3 and R = 10 together with the func- 
tion (1 — r / R)e~ 1 ' 3r . The maxima are chosen to coincide 
with the maxima of the B-spline basis functions. Dashed 
lines: B-spline basis set of order k — 3 with knot sequence 
{0,0,0,2,4,6,8,10,10,10}. 



and r max . We still have two unknowns at this point. We 
make the choice to choose n and determine a for a given 
n from 



+ n(r n 



R) 



' max \i max 



R) 



(5) 



We now construct the BSTO basis by choosing a value of 
r max and n for each function. In practice we spread out 
the values for r max evenly between and R. The value 
of n determines the diffuseness of the BSTO functions. 

Our choice of the BSTO parameters is best illustrated 
with an example. In Fig. [1] we show a typical B-spline ba- 
sis set with 7 functions of order 3 (dashed lines). B-spline 
basis functions have already been used successfully to ob- 
tain continuum properties of atoms and molecules (23[. 
In the same figure we plot a BSTO basis set, where the 
r max are chose to coincide with the B-spline maxima. 
The similarity between the BSTO and B-spline basis set 
is clear. An obvious difference is that the B-spline func- 
tions are only defined between the knots and are zero 
outside, while the BSTO functions are nonzero in the 
entire range from to R. Based on the similarities be- 
tween our BSTO basis set and the B-spline functions we 
expect the BSTO basis set to be successful in describing 
continuum states just as the B-spline basis. 



C. Properties from TDDFT 

In this section we explain how to obtain bound-state 
and scattering properties from TDDFT [H, [H, when 



a hard wall is imposed. The properties considered are 
the polarizability and electron-scattering phase shift. 

The TDDFT equations are usually solved with the lin- 
ear response regime. The key quantity in this regime is 
the frequency dependent first order induced density 



(6) 



The induced dipole moment, which is directly related to 
the polarizability, can be obtained from the linear re- 
sponse density as 



8ji{uj) = — J r5p(r,cj)dr 
and is related to the polarizability by 



(7) 



(8) 



To obtain the induced density we need to know the in- 
duced Kohn-Sham (KS) potential, 5v s (r, a;), and the KS 
response function Xs- The induced KS potential con- 
sists of the induced Hartree potential and an exchange- 
correlation (xc) potential. The induced xc-potential is 
obtained from the xc-kernel, /xcfr,^,^), for which the 
simplest choice is the adiabatic local density approxima- 
tion (ALDA), 



^j2 p unif 

5^{v-v')5{t-t') xc 



dp a dp y 



(9) 



where a and r are spin indices. The density-density re- 
sponse function, Xs > is obtained from 



Xs(r,r',u;) = ^(fi - fa) 



0*(r)0 Q (r)^(rQ^(rQ 
(Si -e a )+^ + irj 



. (10) 



In this expression, i and a indicate occupied and virtual 
states respectively, the fa are the occupation numbers, 
4>i(r) are the KS orbitals, and the corresponding KS 
orbital energies. The response function is completely de- 
termined by ground-state DFT orbitals and orbital en- 
ergies. All these equations are valid for finite systems 
and can still be applied when the density and orbitals 
are confined by a hard wall potential. 

To obtain electron scattering phase shifts we bind an 
extra electron to the A/"-electron system, do a DFT cal- 
culation on this (TV + l)-electron system, and subse- 
quently excite this electron into the continuum using 
TDDFT [28], or, equivalently, excite the electron into 
a discretized continuum state. We obtain the transition 
energies by solving the linear response equations as de- 
rived by Casida [29|. These equation can be written in 
the form of a Hermitian eigenvalue equation: 



n S/T F S/T = 2 p 



S/T 



(ii) 



where the uji are the excitation energies and S and T 
stand for singlet and triplet respectively. The ft matrices 
are given by, 

ttfajb = SabSijiSa-Si) 2 (12) 

+ 2 Vi^" (Kll jh + K} a \ jb ) y^I] 



lL ia,jb 



(13) 



TT 

iajb 



K n 



) 



where the Si and Sj are occupied KS orbitals and e a and 
Sb are unoccupied KS orbitals. The coupling matrix K 
is given by, 



<S>i<r{ v )<t>a<j(*) 



1 



Jr-r'] 

+f^(r,r',u;)}4> b r(r')4> j r(r')drdr', 



(14) 



where is the unknown exchange-correlation (xc) ker- 
nel, for which we use the ALDA kernel. Once the singlet 
and triplet transition energies are obtained, the phase 
shifts S n i for a short-ranged potential can be obtained 
from [5[ 



tan S n i 



Ji(knR) 
ni(k n RY 



(15) 



where ji and n\ are the spherical Bessel and Neumann 
functions, R the hard wall radius, and k n = ^2{uj n — 7), 
where / is the ionization energy. For a long-ranged po- 
tential we have Q 



tan S n i = 



Fi(v,k n R) 
Gi( V ,k n R) : 



(16) 



where Fi and G\ are the regular and irregular Coulomb 
functions and r] is a Coulomb parameter. 



III. COMPUTATIONAL DETAILS 

A. Confined atoms 

In this paper we use the ADF program package [50, HH 
for all confined atom calculations. We use the standard 
ADF even-tempered ET-pVQZ basis set to generate con- 
fined BSTO functions. The basis set used to fit the den- 
sity [24] consists of a large number of ordinary STOs, 
instead of BSTOs. Because of the large number of fit 
functions included, the confined density is in general well 
described by ordinary STOs. When the confinement ra- 
dius becomes very small, the diffuse fit functions can lead 
to numerical instabilities and the most diffuse functions 
are either removed or replaced by more contracted func- 
tions to accurately describe the density of the confined 
atoms. 



The ground-state xc-potential we use for the confined 
atom calculations is the statistical average of orbitals po- 
tential (SAOP). The SAOP potential is asymptotically 
correct and has proven to give good HOMO orbital en- 
ergies and response properties [32[. We note that for 
confined atoms it is not possible to use the near exact 
He [33| and Ne [34| potentials since these potentials cor- 
responds to the density of the free, unconfined, atoms. 
The accurate potentials have a long ranged asymptotic 
tail, while the exact potential corresponding to the con- 
fined atoms is zero at the wall. The exact potentials are 
therefore only asymptotically correct when we do not im- 
pose a wall. 

To obtain the TDDFT polarizabilities we use the 
ALDA kernel. 



B. Scattering phase shifts 

For the BSTO phase shift calculations we use the ADF 
program package [HO, [HI- A BSTO adapted standard 
ADF ET-pVQZ basis is augmented with BSTO functions 
that span the continuum as explained in section Hi Bl As 
we mentioned, the choice of n and r max is, in principle, 
completely free. We let the user choose a value n sta rt > 2 
and then determine the values of all M basis functions 
rii (i = 1, 2, . . . , M) as, 



^start + (i ~ 1) 



50 — n sta rt 

M 



(17) 



where the brackets, [1, indicate that we round to the 
nearest integer of the expression on the right. The choice 
of the riM value to always be smaller than 50 is based 
purely on the limits of the ADF program. The choice of 
letting rii increase steadily, generates functions of more 
even width and makes sure that the value of a does not 
become very small, which is numerically desirable. The 
"width" of the BSTO functions can be modified by the 
choice of n sta rt- Linear dependency problem for larger 
number of functions M can often be solved by choosing 
a larger value of n sta rt- 

We spread the maxima equally in the space between 
2 bohr and R. The reason to start at 2 bohr is that 
the ET-pVQZ basis already spans the region close to the 
nucleus completely. Adding more basis functions in that 
region only leads to linear dependencies in the basis set. 
For the atoms studied in this paper, H, He, and Ne, 2 
bohr turns out to be a suitable value. 

Since we place the wall outside of the short-range part 
of the potential, we can, in principle, use the available 
near exact ground state potentials for the (A^+l)-electron 
systems, H~ and He. The error we make in the potential 
at the radius R is exponentially small if we choose R 
large. In practice we only use the near exact potential 
for He [33|]. In case of H~ (e-H~ scattering) the near 
exact potential is not available at enough grid points to 
obtain accurate results. In Ref. [4] it was shown that 
for this system the EXX potential is very close to the 
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FIG. 2: HOMO orbital energies of He and Ne as a function 
of the wall radii, compared with HF results from Ref. Q- 
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FIG. 3: HOMO orbital energies of He and Ne as a function 
of the wall radii, compared with HF results from Ref. 0. 



exact potential. Therefore we use the EXX potential as 
input for our calculations. We apply the EXX potential 
as a fixed input potential, even though the ADF code 
has the option to perform Hartree-Fock exchange-only 
and hybrid potential calculations. The reason is that 
our BSTO method is, as yet, not implemented for the 
case of these orbital dependent potentials. The orbital 
dependence of the potentials is precisely the reason that 
implementation of the BSTO basis set for these potentials 
is considerably involved and left for future applications. 

We compare the BSTO results with results obtained 
with a well-established fully numerical spherical DFT 
code, which includes the optimized effective potential 
method (OEP) and is supplemented by the option to 
insert a hard- wall at a distance R from the origin Q. 
This program is basis- set independent, works with a ra- 
dial grid, and both the energies and the potentials are 
optimized in a self-consistent way. The TDDFT exci- 
tation energies are calculated with a supplemental code 
that explicitly solves the radial TDDFT equations [4(. 

We use the ALDA kernel for all TDDFT calculations. 



IV. RESULTS 
A. Confined atoms: He and Ne 

In this section we show results for the HOMO energy 
and the static polarizability as a function of the box ra- 
dius for two many-electron atoms, He and Ne. In Fig. [2] 
we show the change in HOMO energy with the wall radius 
for both He and Ne. The HOMO KS orbital decays expo- 
nentially and the wall will therefore have some influence 



up to oo and the free atom value is approached asymp- 
totically. From the figure it is clear that the HOMO 
energies are already very close to the corresponding free 
atom values around 4 bohr. We obtain 95% of the free 
atom He HOMO energy at a radius of 3.06 bohr. For Ne, 
the radius where 95% of the HOMO energy is obtained 
is at 3.50 bohr. 

In the region where the HOMO energies are practically 
leveled off, i.e. after 4 bohr, the difference between the 
HF and SAOP values is almost constant. For He we find 
a free atom value of -0.8165 a.u., the HF reference value 
is -0.9174. For Ne our SAOP free atom energy is -0.7585 
a.u. and the reference HF value is -0.8511 a.u.. The 
differences are a result of the difference between the HF 
potential and our SAOP potential. The free atom value 
of He is equal to the HF value when the exact xc-potential 
is used, but as mentioned in the previous section, this 
potential is not valid when the atom is enclosed in a box 
of finite radius. We can try to correct for the shift by 
vertically shifting the SAOP curves with the free atom 
energy difference. We show the result in Fig. [3j The 
shifted He and Ne results are almost on top of the HF 
results. 

The ionization radius is defined as the radius where 
the HOMO level of the confined atom jumps from a neg- 
ative to a positive energy when the box becomes smaller. 
The "ionized" electrons are still confined to the box, and 
therefore not explicitly ionized as such. The atomic ion- 
ization radii are 1.50 bohr for He and 1.89 bohr for Ne, 
these values correspond well with the 1.41 bohr and 1.83 
bohr of Ref. Q. In Ref. [H| a DFT result for the ion- 
ization radius of Ne was obtained with the PW potential 
using total energy differences, they obtained a value of 
1.83 bohr. Our SAOP results are close to the Hartree- 
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FIG. 4: Static dipole polarizability of He and Ne as a function 
of the wall radii. The dotted horizontal lines indicate the 
respective free atom values. 



Fock (HF) values of Ref. 0. When we use the shifted 
SAOP curves to recalculate the ionization radii we obtain 
1.45 bohr for He and 1.84 for Ne, very close to the HF 
values. 

In this work we present, to our knowledge, the 
first TDDFT polarizabilities for confined many electron 
atoms, and the first results including correlation. All 
previous works on many electron atoms used the HF ap- 
proximation and obtained the polarizabilities using fi- 
nite field methods. We use the SAOP potential in the 
ground state and the ALDA kernel for the response cal- 
culation of the polarizability. We show our results for He 
and Ne in Fig. [H For small wall radii we see that the 
polarizability approaches zero. This is consistent with 
the findings that for one-electron systems, the results 
of which are expected to at least qualitatively apply to 
many-electron systems, the polarizability increases from 
R = as ~ R 4 [36[. For large radii the polarizability 
asymptotically approaches the free atom value from be- 
low as expected. 

The free atom values we predict are 1.425 for He and 
2.565 a|j for Ne, which correspond well with the experi- 
mental values of 1.41 clq for He and 2.57 clq for Ne [37( ] . 
We obtain an overall s-shaped curve, in agreement with 
the s-shaped curves found in Ref. [38|, [39[ for several 
many electron atoms within the Hartree- Fock- Slater ap- 
proximation. It is clear that the asymptotic approach of 
the Ne polarizability towards the infinite energy values 
is much slower than that of the He atom. As can be 
seen from Eq. [10l the polarizability depends on all oc- 
cupied and unoccupied states, including the continuum. 
The wall always has an influence on the higher Rydberg 
states, since the potential of the neutral atoms is long- 
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FIG. 5: e-He + s- and p-wave scattering phase shifts (bare 
KS, singlet, and triplet) for BSTO basis sets with different 
wall radii. The bottom three curves correspond to p-wave 
scattering and the upper three curves to s-wave scattering. 
We compare the BSTO results to the phase shifts for a large 
radius of 100 bohr. Reference literature results are also in- 
cluded 53,51. 



ranged. The (pseudo) continuum states in the box con- 
verge when the short-ranged part of the potential, that 
part that differs from — 1/r, has converged. For He this 
point is reached well before 10 bohr, while for Ne this 
point is reached much later. This explains why the con- 
vergence of the Ne polarizability is so much slower than 
the He polarizability. 

We expect the results for He and Ne to be representa- 
tive for other closed shell atoms. For both the HOMO 
energies and the polarizabilities we reproduce the trends 
observed in other works. The BSTO basis set works well 
for these systems. 



B. Electron scattering: e-He + and e-H 

We show results for the phase shifts corresponding to 
e-He + and e-H scattering in Figs. [5] and [6l We compare 
our BSTO results with grid calculations corresponding 
to a wall radius of R = 100 bohr. The results of the 
different methods should coincide in case of a large basis 
set and a large grid. 

We first need to decide where to place the wall. We 
want to place the wall as close as possible to the nucleus 
to limit the necessary basis set size, but it needs to be in 
a region where the density and the occupied orbitals are 
negligibly small In case of He, the density has a range 
of approximately 2 bohr and the HOMO a range of ap- 
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FIG. 6: e-H s-wave scattering phase shifts (bare KS, singlet, 
and triplet) for BSTO basis sets with different wall radii. We 
compare the BSTO results to the phase shifts for a large ra- 
dius of 100 bohr. Reference literature results are also in- 
cluded m. 



proximately 6 bohr. We place the wall ati? = 10,ll,12 
bohr for this system. In case of H~ , the density also has 
a range of approximately 2 bohr, but the HOMO orbital 
is very diffuse (its energy is only -1.257 eV) and it has a 
range of approximately 20 bohr. Therefore we place the 
wall at R = 20, 21, 22 bohr. 

In Fig. [5] we see that results obtained with the BSTO 
basis set are indeed equal to the grid code results for 
e-He + with I = 0,1. Because we only have access to 
relatively small radii using the BSTO basis set, the points 
are far apart (the energies are widely spaced). To obtain 
the full curve, many different wall radii need too be used. 
Since the calculations are computationally very fast, this 
is in general not a problem. We see that even with only 
three different radii the curve is already well represented. 
On average 12 extra s- and p-functions are necessary on 
top of the standard ADF ET-pVQZ basis set to obtain 
these results. 

In Fig. [6] we show the BSTO results together with large 
radius results for e-H scattering with 1 = 0. Since the wall 
distance is larger, we sample more points for each range. 
Again the entire curve is well reproduced with only 3 wall 
radii. Only 15 extra s-functions were necessary on top 
of the standard ADF ET-pVQZ basis set to describe the 
range of 0-22 bohr. 

Overall, we observe that the phase shifts are well re- 



produced by the BSTO basis set method. 

V. CONCLUSIONS 

In this paper we placed an atom inside an infinite 
spherical potential by using a modified STO basis set 
multiplied by a cutoff function. We used (TD)DFT for 
the calculation of the HOMO energy and polarizability in 
case of confined atoms and for the calculation of electron- 
scattering phase shifts. 

By placing the potential wall at different distances 
close to the atomic nucleus we studied the evolution of 
the HOMO energy and polarizability of the confined He 
and Ne atoms. The HOMO energy and polarizability ob- 
tained with (TD)DFT follow the same trends as available 
HF results @, [38|, [3§] . 

By placing the wall at a large distance from the 
atomic nucleus, we discretized the continuum and ob- 
tained electron-scattering phase shifts for e-He + and e-H 
scattering. We showed that the same accuracy can be 
obtained with the basis set method as with a recently es- 
tablished grid method The systems studied in this 
paper do not exhibit any low-energy resonances. When 
resonances are present the basis set method can be used 
in conjunction with the stabilization method [43[ to find 
their locations. 

This research paves the way to study molecules in cav- 
ities and electron-molecule scattering using TDDFT. Go- 
ing from the spherically symmetric system studied in this 
paper to molecules is not trivial. First we need to de- 
scribe atoms that are off-center in a spherical cavity, this 
makes the evaluation of 1-center integrals more compli- 
cated. Such integrals have for example been evaluated for 
an off-center He atom in a spherical cavity using GTOs 
in Ref. [45[. In case of molecules we also need to take 2- 
center integrals into account, a way to treat this problem 
using floating GTOs is described in Ref. [46]. Finally, 
the theory of electron-molecule scattering using TDDFT 
needs further development. 

In conclusion, we showed that a unified basis set 
method can be applied to study both (time-dependent) 
properties of confined atoms and low-energy electron- 
atom scattering. 
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